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ABSTRACT 



Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high- 
z AGN hosts. Combining PACS 100 and 160yum observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% 
of X-ray AGN individually at > 3it. The host far-infrared luminosity of AGN with L2-iokev ~ 10*^^ erg s"' increases with redshift 
by an order of magnitude from z=0 to z~l. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for 
L2-iokcV S lO^'^ergs"' AGNatz>l. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample 
which is dominated by L2_iokcV < 10''"' ergs"' AGN. In conjunction with properties of local and luminous high-z AGN, we interpret 
these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star 
formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects 
an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive 
galaxies and shows littie dependence on AGN luminosity. The levei of this secular, non-merger driven star formation increasingly 
dominates over the correlation at increasing redshift. 

Key words. Galaxies: active - Infrared: galaxies 



1. Introduction 

Measuring the star formation rate of the host galaxy is im- 
portant for studying the co-evolution of active galactic nuclei 
(AGN) and their hosts. It is often difficult because the AGN 
can outshine the host at many wavelengths. However, the rest 
frame far-infrared/submm emission appears dominated by the 
host for AGN with vLy(60jum)«0.1. . .0.2 LboI.agn and higher 
(e.g. Netzer et al. 120071 an d introduction to Lutz et al. 120101 
see also Wang et al. 120081 for luminous high-z QSOs) and has 
been used as a host star formation rate diagnostic of high-z AGN 
(e.g. Serjeant & Hatziminaoglou 2009, MuUaney et al. 2010 and 
Lutz et al. l2010l) . Herschel can detect much lower star formation 
rates than previous far-infrared and submm studies. It further im- 
proves such work by measuring the rest frame far-infrared SED 
peak without extrapolation from longer wavelengths. We here 
present a first Herschel study of rest frame far-infrared emission 
and host star formation in X-ray selected AGN in the GOODS-N 
field. Throughout the paper, we adopt an Q,„ = 0.3, Oa - 0.7 
and //() - 70 kms"' Mpc"' cosmology. 

2. Data 

We use the v2.2 100 and 160yum ma ps of t he GOODS-N field 
obtained with PACS (Poglitsch et al. 120 lot on board Herschel 
(Pilbratt et al. 2010) as a science demonstration project for the 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-fed Principal Investigator consortia and with impor- 
tant participation from NASA. 



PACS Evolutionary Probe (PEP) guaranteed time survey. We use 
the PACS catalog extracted with IRAC/MIPS24//m based posi- 
tion priors by B. Magnelli, to a 3cr depth of ~3.0 mJy and ~5.7 
mJy at 100 and 160jum, respectively. Since there are Chandra 
sources that are not detected at 24fim, we have verified by in- 
dividual comparison to a blind catalog that we are not missing 
PACS detections on such sources when using the 24//m prior 
catalog. For samples of objects undetected by PACS we stack, 
using the X-ray positions, into the residual maps obtained af- 
ter subtracting the sources in the prior catalog from the original 
maps. See also the appendix to Berta et al. (2010,) for a descrip- 
tion of the data and catalogs. 

The CDFN 2 Msec Chandra X-ray catalog of Alexander et 
al. (2003) provides the basis for our AGN selection. 328 of its 
503 X-ray sources lie within the part of the PACS map, which 
has at least half of the coverage of the central homogeneously 
covered region; we restrict ourselves to this subset. Reaching 
a: 2.5 X lO^'^ergcm-2 g-i jjj jjje 0.5-2keV band, CDFN X-ray 
detections include a noticeable number of star forming galax- 
ies without significant AGN contributions, in addition to bona- 
fide AGN. Bauer et al. (2004) have outlined for this very sample 
a combination of criteria to distinguish X-ray AGN from star 
formation dominated objects, based on X-ray luminosity. X-ray 
obscuring column or hardness, optical spectroscopic classifica- 
tions, and X-ray/optical flux ratio. We use an updated version 
of their classifications for the GOODS-N X-ray sources, which 
separates the 328 sources with good PACS coverage into 224 
X-ray AGN, 67 X-ray detected star forming galaxies, and 37 
other targets (stars or unidentified). We use both spectroscopic 
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Table 1. PACS detections of GOODS-N X-ray sources 



Class 


Sources 




> 3(T detections at 






lOOfim 


160yum 


100 


or 160yum 


All 


328 


76 


81 




91 


Galaxies 


67 


37 


37 




40 


AGN 


224 


35 


39 




46 


Other 


37 


4 


5 




5 



(Barger et al. l2008l 57% of our X-ray AGN) and photometric 
redshifts (F. Bauer, mostly based on Barger et al. 2003). X-ray 
spectral fitting (F. Bauer et al., in prep.) provides intrinsic X-ray 
luminosities L2-i()kev and obscuring columns Nh for the X-ray 
sources. The fits use absorbed powerlaws with a fixed Galactic 
column of 1 .6x10^" cm"^ and a variable obscuring column at the 
redshift of the source. The photon index is allowed to vary for 
sources above 150 net counts in the 0.5-8keV band, and is fixed 
to 1.85 below. 
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Fig. 1. Redshift distribution of the 224 GOODS-N X-ray AGN 
in the region with good PACS coverage. 46/224 are individually 
detected in at least one of the PACS bands. 



3. Results 

The far-infrared detection rate of the 328 X-ray sources above 
the 3cr level in at least one of the two PACS bands is 28% 
(Table [U. For the three subgroups, it is 60% (galaxies), 21% 
(AGN) and 14% (other). The large detection rate of 'galaxies' is 
consistent with their X-ray emission arising from star formation 
related processes. A detailed discussion will be provided else- 
where. In the following, we focus on the 224 X-ray AGN only. 

We use both the individual detections and stacks of non- 
detections derived via the stacking library of Bethermin et al. 
(12010). The 178 AGN that are individually undetected in both 
PACS bands are statistically well detected in the stacks at 
0.79+0.08mJy (100//m) and 1.46+0.16mJy (160yum). These de- 
tection errors are derived from multiple stacks at random posi- 
tions over the same coverage region of the PACS residual maps, 
and using the number of targets as the sample of interest. 

To derive far-infrared (FIR) luminosities with minimal as- 
sumptions about SED shapewe compute the rest frame 60jum 
luminosity vLv(60;um) using the detection wavelength closer to 
rest 60/jm or log-linearly interpolating for detections in both 
bands at 0.67<z<1.67. We treat stacked 100/im and 160/im 
fluxes equivalently to derive the mean rest frame 60/im lumi- 
nosities of the individually undetected sources. Given the strong 
K-corrections for PACS fluxes over the redshift range of X-ray 
AGN, this requires restricted redshift ranges for the stacks, we 
limit to bins with Az/(1 + z) < 0.4 and adopt the median z of 
the particular sample. Luminosities for the combined sample of 
detections and nondetections are obtained as averages weighted 
by the number of sources. 

Figure[T]shows the redshift distribution of all X-ray AGN in 
our sample. As expected, PACS photometry is most efficient in 
detecting z<2.5 AGN hosts. We also focus on z<2.5 in order to 
probe far-infrared rest frame > 45/vm wavelengths, beyond the 
mid-infrared that is more easily AGN dominated. Table |2] and 
Fig. |2] show 60//m luminosities as a function of redshift, sepa- 
rately for the PACS detections, stacks of nondetections and the 
combined sample. Host 60//m luminosities increase with red- 
shift. An increase is seen in detections, stacked nondetections, 
and in the averages for the combined sample. This trend cannot 
be simply due to the increase of FIR detection limit with redshift, 
which would leave the average luminosities from the combina- 
tion of detections and nondetections unchanged. 
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Fig. 2. Far-infrared luminosities of GOODS-N X-ray AGN as a 
function of redshift. Symbols labelled 'AH AGN', with the red- 
shift range indicated, average the detections and nondetections. 
Their uncertainty is derived from bootstrapping into the com- 
bined sample. Number of detections and total number of sources 
in each bin are indicated. The insert shows the stack for the ind- 
vidually undetected 0.8<z< 1.4 sources. 



Mean fluxes for the PACS > 3cr detections and stacked mean 
for the nondetections differ by about an order of magnitude in 
each of the three redshift bins. Such a large difference must re- 
flect a wide intrinsic distribution of far-infrared luminosities. For 
the simplifying assumption of a log-normal distribution of far- 
infrared luminosities, the ratio of mean detections to the stacked 
mean as well as the ~20% detection rates can be reproduced with 
an intrinsic dispersion of about 0.5 dex for the distribution of far- 
infrared luminosities in each bin. The detailed shape of this dis- 
tribution is not constrained but must be wide. The mean 160/vm 
flux for all the 1.4<z<2.5 AGN is 3mJy (in agreement with ten- 
tative MIPS values for ECDFS AGN by Papovich et al. 2007}. 
Given the wide distribution of FIR luminosities, the typical (me- 
dian) flux must be lower, by a factor ~2 for the log-normal dis- 
tribution. A considerable variety exists not only in far-infrared 
luminosity but also in mid-to far-infrared SED shape (Fig. O, 
from SEDs quite similar to those of star forming galaxies to ones 
that are clearly AGN dominated over a wide wavelength range. 



Table 2. Mean FIR luminosities of different AGN groups 
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z 


L2-10kcV 


Noe./N 


Detected 


Stack 


All 






erg s"' 




vL^(60yum rest), 


10*" erg S-' 


0.2- 


-0.8 


all 


12/35 


3.3 


0.32 


1.33±0.31 


0.8- 


-1.4 


all 


22/83 


11.3 


1.43 


4.04±0.62 


1.4- 


-2.5 


all 


10/76 


61.1 


5.27 


12.61+3.12 


0.2- 


-0.8 


<10^^ 


4/16 


1.9 


0.09 


0.55±0.23 


0.8- 


-1.4 


< lO''^ 


1/9 


6.9 


1.40 


2.01±0.91 


1.4- 


-2.5 


<10'*2 


0/1 




3.92 




0.2- 


-0.8 


1042-4J 


2/12 


4.4 


0.60 


1.23±0.51 


0.8- 


-1.4 


1042-43 


7/37 


12.3 


1.56 


3.60±0.73 


1.4- 


-2.5 


1042-43 


6/18 


63.4 


4.35 


24.04±8.45 


0.2- 


-0.8 


1043-44 


6/7 


3.8 


1.91 


3.54±0.89 


0.8- 


-1.4 


1043-44 


13/32 


10.5 


1.61 


5.22±1.11 


1.4- 


-2.5 


1043-44 


3/51 


69.3 


5.30 


9.06+2.80 


0.2- 


-0.8 


>10^* 


0/0 








0.8- 


-1.4 


>10'"' 


1/5 


18.2 


2.32 


5.50±3.56 


1.4- 


-2.5 


>10'*^ 


1/6 


22.2 


8.58 


10.85±4.77 
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Nh 


NDe./N 


Detected 


Stack 


All 






cm ^ 




vL,,(60yum rest), 


lO-*-* erg S-' 


0.8- 


-1.4 


<10^^ 


2/10 


13.7 


1.12 


4.15±1.25 


0.8- 


-1.4 


1022-23 


6/25 


11.6 


1.47 


3.89±0.91 


0.8- 


-1.4 


1023-24 


10/28 


9.6 


1.76 


4.57±0.89 


0.8- 


-1.4 


> 10-^ 


0/5 




1.67 


1.67+0.88 



Note: Average FIR luminosities are given separately for the individual 

> 3(T detections, the stack of the nondetections, and the 

number- weighted mean for all sources. Errors for luminosities are 

standard deviations from bootstrapping and are dominated by the 

dispersion of the AGN population rather than measurement error 



















X 








t 


- 




10-" 










X 










* 


); 


: 


F 




























u 










¥ 


y^ ^ m 


■ 












- 




10"" 


r 




M 


X 
■ 


■ 


■ 

■ 
■ 


« 










- 




























^ 








>« 




■ 


















- 






■ A03 


ID: 256 


L2-io«.v=2.8x10" 


erq 


s"' 




10-" 


r 


* 


■ 






X A03 


ID: 390 


'-2-10 keV-9-4x 


lO*' 


erg 


s~^ 




- 




























, 




. 






, , 


, , 




1 - 



1.0 10.0 

Rest-frame wavelength [fJ.m] 



Fig. 3. Examples of z~1.2 AGN illustrating the range of optical 
to far-infrared SEDs. ID 390 (Alexander et al. |20U3] ) is domi- 
nated by AGN continuum over the optical to mid-infrared, while 
ID 266 which hosts an high obscuring column AGN is dom- 
inated by an optical/NIR stellar bump and a pronounced far- 
infrared peak. 



The increase of host far-infrared emission with redshift in 
Fig.|2]could still be influenced by the increase of AGN luminos- 
ity with redshift that is inherent to this X-ray selected sample. To 
break possible degeneracies, we have binned our sample in 3x4 
bins defined by redshift (z=0.2-0.8, 0.8-1.4, 1.4-2.5) and by in- 
trinsic hard X-ray luminosity (L2-iokeV< 10"*^ergs"', lO'^^-lO'*'', 



10«_io44, > 10^4). Table|2]Hsts FIR luminosities for these bins, 
for the figures in the following discussion we omit bins with less 
than 5 objects and correspondingly large errors. Errors on the 
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Fig. 4. Left; Far-infrared luminosity as a function of redshift, 
for different bins in intrinsic rest frame 2-lOkeV X-ray lumi- 
nosity. Values reflect the mean of detections and nondetections, 
and errors are based on bootstrapping into the respective sam- 
ple. Right: FIR luminosity as a function of intrinsic hard X-ray 
luminosity, for different redshift bins 
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Fig. 5. Far-infrared luminosity as a function of X-ray obscuring 
column. Sample forredshifts 0.8<z<1.4. Low obscuring column 
objects have been placed at 2 x 10^"cm"^. Symbols are as in 

Fig.ia 



average FIR luminosity are dominated by the variations in the 
underlying population rather than by photometric error. For that 
reason, errors in Table |2] and the figures are standard deviations 
derived from bootstrap estimates for the respective subsamples. 
Fig |4] (left) shows that the increase of host FIR luminosity 
with redshift is clearly preserved when considering AGN lumi- 
nosity bins separately. Focusing on the far-infrared luminosities 
of L2-iokev = l 0^*^-1 0'*^ AGN, for which there are more than 10 
objects in each redshift bin, there is no overlap in the 99% confi- 
dence intervals of FIR luminosity comparing the z=0.2-0.8 and 
the z=1.4-2.5 redshift range (< 2.8 vs. > 6.1 x lO'^'^ergs"'). 
These 99% confidence intervals were directly derived by boot- 
strapping and thus include non-gaussianity of distributions of 
individual source properties or of errors. The higher the AGN 
luminosity bin we consider, the higher is the FIR luminosity at 
low z but then its increase with redshift less steep. Conversely, 
when studying FIR luminosity as a function of AGN luminosity 
(Fig. |4] right) there is an increase with AGN luminosity in the 
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results in that context. Here we have converted from hard X-ray 
to AGN bolometric luminosity using eq. (5) of Maiolino et al. 
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Fig. 6. Star forming (=far-infrared) luminosity vs. AGN lumi- 
nosity for the GOODS-N AGN and a local reference sample 
of extremely hard X-ray selected BAT AGN. The dotted col- 
ored lines indicate schematically how the observations are ex- 
plained by the combination of a diagonal 'evolutionary connec- 
tion' trend with a general increase of host star formation with 
redshift in hosts of moderate luminosity AGN, similar to that 
for the general galaxy population. The dashed line is the relation 
imphed by Netzer et al. (2009) . 

lowest redshift bin that flattens at higher redshift, with no signif- 
icant trend left at z>1.4 . 

In the luminosity range covered by our sample, we do not 
find a significant trend of FIR luminosity with X-ray obscur- 
ing column (Fig. |5]). Such a trend would be expected in merger 
evolutionary scenarios that are invoking a sequence starburst - 
obscured AG N - un obscured AGN (e.g. Sanders et al. 119881 
Hopkins etal. l2006] l. 

4. Evolution of the relation between AGN luminosity 
and host star formation 

Due to the excellent sensitivity of HerschelPACS to rest frame 
far-infrared emission in the hosts of z=0.2-2.5 X-ray selected 
AGN, we have been able to break the redshift/luminosity degen- 
eracy that aff'ected previous submm-based studies of AGN host 
star formation (e.g. Lutz et al. 1201 Oi l. We avoid having to make 
SED assumptions that were necessary for submm-based stud- 
ies. Compared to the Spitzer 70//m-based study of Mullaney et 
al. (2010), we are probing the rest frame far-infrared out to z>l 
where 70/im data already enter the rest frame mid-IR with less 
favourable contrast between host emission and AGN heated dust. 
Results of these submm as well as 70;um based studies agree well 
with our findings within the mentioned constraints. To study evo- 
lution over a yet wider redshift range, we supplement in the fol- 
lowing our GOODS-N sample with local (z< 0.3) Swift B AT 14 - 
ISOkeV extremely hard X-ray detected AGN (Cusumano l2009l l. 
Here, far-infrared information is taken from the IRAS database, 
in a way consistent with the treatment of the GOODS-N AGN 
(see also Sect. 3.3 of Lutz et al. 2010). 

A variety of studies have established a local correlation 
between AGN luminosity and star forming luminosity (e.g. 
Rowan-Robinson 1995; Netzer et al. ' 20071 |2TO I I. at least in the 
luminosity regime of bright Seyferts or Quasars. High redshift 
QSOs appear to extend this relation to yet higher luminosities 
(e.g. Lutz et al. l2008l and references therein). Figure |6]places our 



(I2007I I and a ratio 7 between bolometric and 5100A luminosity. 
Local BAT AGN follow the coiTelation of AGN and star forma- 
tion luminosity, with the exception of a flattening at low AGN lu- 
minosities. Figure |6] shows that high luminosity AGN stay close 
to the correlation for all redshifts covered by our sample, but the 
host star formation rates of lower luminosity AGN rise by about 
an order of magnitude from z=0 to z=l and by about 1.5 or- 
ders of magnitude from z=0 to z=2. Such a behavour can be ex- 
plained by the combination of two paths of AGN growth. On one 
path, AGN growth and host star formation are tightly coupled 
by an evolutionary mechanism, likely merging. This will result 
in the diagonal correlation line in Fig. |6] The other path reflects 
a secular evolution with no close coupling of AGN luminosity 
and galaxy-integrated host star formation rate. Connections be- 
tween AGN and star formation phenomena on smaller spatial 
scales or diff'erent timescales might however be present. Here, 
star formation levels will be similar to those that are pervasive 
in massive galaxies at a given redshift. This coiTesponds to the 
low-luminosity flattening of the relation in Fig. |6]at a level in- 
creasing with redshift. The detailed slope of this flatter part is 
still poorly constrained at current statistics. Bouche et al. (1201 Oi l 
have parametrized on the basis of a variety of observational stud- 
ies star formation rates of star forming galaxies as a function of 
galaxy mass and redshift. Their equation (1) corresponds to an 
increase in SFR from redshift to the centers of our three red- 
shift bins by factors 3.0, 7.4, and 19 respectively, consistent with 
the location of AGN on the flat 'secular' path in Fig. |6] As dis- 
cussed in Lutz et al. ( 1201 Oi l, such a two-path scenario for AGN 
is also consistent with a variety of other properties of AGN hosts 
at high z. 

A limitation to the current study is the small ~ 11' x 16' 
field which restricts the number of L2-ioiieV > lO^'^erg s"' AGN. 
Upcoming Herschel observations will provide results for larger 
samples of such objects, where we will be able to investigate the 
relation above these AGN luminosities, i.e. on the 'merger path'. 
We might expect an upturn in host far-infrared luminosity which 
is suggested at zwl for more luminous AGN (Lutz et 1201 Oi l, or 
the possible relation of host star formation to AGN obscuration 
(e.g. Page etal. l200TT l. 
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